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Abstract 

It is shown that the two qubit CNOT (controlled NOT) gate 
can also be realised using q-deformed angular momentum states con- 
structed via the Jordan-Schwinger mechanism. Thus all the three gates 
necessary for universality i.e.Hadamard, Phase Shift and the two qubit 
CNOT gate are realisable with q-deformed oscillators. 
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1. Introduction 

Recently^ it has been shown that the single qubit quantum logic gates, 
viz. ,the Hadamard and Phase Shift gates can also be realised with two q- 
deformed oscillators where q is the deformation parameter of a quantum Lie 
algebra^. q-Deformed oscillators here mean that the Lie algebra satisfied by 
creation {a'^) and annihilation (a) operators of a bosonic oscillator, viz. aa^ — 
a'^a = 1 is modified into a^aj — gaja^ = where N is the number operator 
and q — ; 0<s<l. Using such deformed oscillators an alternative 
formalism for quantum computation can be set up^. The advantage of this 
over the conventional formalism (which is obtained for g ^ 1) is the presence 
of an arbitrary function which may be exploited for experimental purposes. 

However, the formalism will be more meaningful if the realisation with q- 
deformed qubits is possible for all the gates required for universality. A set of 
gates is said to be universal for quantum computation if any unitary operation 
may be approximated to arbitrary accuracy by a quantum circuit involving 
those gates. In the case of standard quantum computation the Hadamard, 
Phase Shift and the CNOT (controlled NOT) gates constitute such a set^. In 
this paper 1 show that the 2-qubit controlled-NOT gate can also be realised 
with q-dcformcd qubits. Thus all the three gates, i.e. Hadamard, Phase Shift 
and CNOT gates can now be obtained with q-deformed qubits. 

The motivation for considering q-deformcd oscillators in quantum com- 
putation comes from the fact that deformed oscillators have been successfully 
used as a tool to understand deviations from an ideal theoretical or exper- 
imental scenario^ for past many years. Bonatsos and Daskoloyamis^ were 
among the firsts to show that the vibration spectra of diatomic molecules 
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gave better fits using deformed oscillators. Parisi"^ studied a d- dimensional 
array of Josephson junctions in a magnetic field and computed the thermo- 
dynamic properties in the high temperature region for o? — > oo. Evalua- 
tion of the high temperature expansion coefficients were done by mapping 
onto the computation of some matrix elements for the q-deformed harmonic 
oscillator. Raychev^ et al calculated the deviations from the nuclear shell 
model using the q-deformed three dimensional harmonic oscillator. Bonatsos, 
Lewis, Ray chev and Terziev^ demonstrated that the three dimensional q- 
deformed harmonic oscillator correctly predicts the first supershell closure in 
alkali clusters without introducing additional parmctcrs.McEWan and Freer^ 
showed that the nuclear orbitals of certain nuclei were commensurate with 
the energy level scheme of the deformed harmonic oscillator and the Nils- 
son model. For these reasons it is meaningful to study whether deformed 
oscillators can be used in the formalism of quantum computation. 

A brief review of relevant facts is given in section 2. In Section 3 the NOT 
gate is realised with q-deformed oscillators. Section 4 gives the realisation of 
the two qubit CNOT gate in terms of q-deformed oscillators. In Section 5 the 
states are discused and Section 6 is the conclusion where a brief elaboration 
of the possibility of quantm error correction using deformed oscillators is 
given. 

2. Brief Review 

Quantum logic gates are basically unitary operators Three gates, the 
single qubit Hadamard and Phase Shift gates and the 2-qubit CNOT gate, 
are sufficient to construct any unitary operation on a single qubit ^.This is the 
universality referred to above. These gates are constructed using the "spin 
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up" and "spin down" states oiSU{2) angular momentum i.e., the two possible 
states of a qubit are usually represented by "spin up" and "spin down" states. 
This is the Jordan-Schwinger construction using two independent harmonic 
oscillators. A similar construction of q-deformed angular momentum states 
can be done using two q-deformed oscillators^^. In Ref.l it was shown that 
the Hadamard and Phase Shift gates can also be realised with q-deformed 
qubits. To achieve this ,the technique of harmonic oscillator realisation ^^'^^ 
of q-oscillators was used. This allows one to set up an alternate quantum 
computation formalism. 

q-Oscillators are described by deformed creation and annihilation opera- 
tors, alj,aq ,respectively. For ordinary oscillators these are a'' and a. q — e*, 
< s < 1 and the deformed oscillators satisfy the following relations : 

a,4-^a, = g-^ ; TVt = TV (la) 

[TV, a,] = -a, ; [TV, aj] = aj ; aja, = [TV]; a,aj = [TV + 1] (16) 

a,/(TV) = /(TV + l)a, ; at/(TV) = /(TV - 1)^ (Ic) 

where the q-number [x] — {q^ — q~^) /{q — q~^) becomes the ordinary number 
X when g — > 1 (i.e.s — > 0). TV is the number operator for the q-deformed 
oscillators and /(TV) is any function of TV. The eigenvalue n of the number 
operator TV denotes the number of bosonic particles. We confine to real q. 
Qg, a J and a, are related as ^-"^ 



ttq = a 



• n I 



\ TV(g-g-i) ' \ TV(g-g-i) 

N = N - {l/s)lnilj2 (26) 



(2a) 
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is the number operator for usual oscillators with eigenvalue h; and ijji , '02 
are arbitrary functions of q only with ■0i_2(?) = 1 for g = 1. 



operator of the deformed oscillator becomes identical to the number oper- 
ator of the standard harmonic oscillator. But oscillator states are usually 
expressed in occupation number basis. So if the number operators are iden- 
tical, there is no way of differentiating between a deformed oscillator state 
and a standard oscillator state. So nothing is gained and we are still in the 
realm of standard quantum computation. But equations (2a, b) are general 
if the ipi{q),i — 1,2 are not all equal to unity. Let ipi — ip2 — i^{q) ■ Now 
N — N — (l/s) ln'4>{q) (equation (26)). This will be reflected in the Jordan- 
Schwinger construction of angular momentum states and the states in the 
two cases will be distinguishable through the function ip{q)- Further details 
are in Ref.[l]. 

We now express a single qubit state in terms of two harmonic oscilator 
states using the Jordan- Schwinger construction. 

(a) States are defined by the total angular momentum j and ^-component 
of angular momentum i.e. m. A particular (j, m) state is created by 
acting the creation operators on the vacuum or ground state which is a direct 
product state of the individual oscillator ground states : 



10) = |0) = |0)i|0)2 is the ground state (j = 0,m = 0), while |0)j,i = 1,2 
are the oscillator ground states, j — {n\ -\- n2)/2 ; m — {rii — n2)/2 where 
ni, n2 are the eigenvalues of the number operators of the two oscillators. 



If all these arbitrary functions are unity, then N = N, i.e. the number 




(3) 
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1^ 1 /o\ l"-2; in\ . \ 1 low l"2; in\ 

- '/^^ = [K)!(l-nO!]V2'°^ > I - (^1 - V2)> = [(^^),(i_^^),]i/2lO) 



(b) A qubit can be either "up" or "down" i.e. there are two possible 
configurations. So the oscillator number operators can take the following 
sets of values only : (rii = l,n2 = 0, "up" state) and (rii = 0, 712 = 1, 
"down" state). Hence {n-i + n2)/2 = 1/2. As j — 1/2 for both qubit states, 
suppress j for simplicity of notation : 

I . _ (aI)^/^+-(4)^/^-" I , |_ _ (aI)V^-(4)V^+-" 

[(l/2 + m)!(l/2-m)!]V2l'^/ ' I [(1/2 + m)!(l/2 - m)!]V2 1*^^ 

(4a) 

Equivalently, in terms of ni, 712 these are 

/ 

(46) 

(c) The basis states are (|1) = \up) state and |0) = \down) state ) 

|1) = |l/2, 1/2) = |l/2) = alio) = al|0)i|0)2 = |i)i|0)2 

|0) = |l/2, -1/2) = I - 1/2) = 4|0) = 4|0)i|0)2 = |0)i|i)2 

Thephysical meaning of the notation is as follows. The |1) (spin "up") state 
is constructed out of two oscillator states where the first oscillator state has 
occupation number 1 while the other has occupation number 0. The |0) ( spin 
" down" ) state corresponds to the first oscillator having occupation number 
and the second oscillator having occupation number 1. So any qubit state 
\x) is : 

k) = (4)^(4)^-16) (5) 

|0) represents one of the two possible qubit states while |0) represents oscilla- 
tor ground state i.e. occupation number 0; |1) represents an oscillator state 
with occupation number 1 etc. 
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Now consider the Hadamard transformation. For a standard qubit, the 
Hadamard gate acts as follows. When one of the basis states is given as an 
input, the output is a superposition of the two basis states, i.e. |0) |0) + 11) 
and |1) |0) — So the Hadamard transformation on a single qubit state 
{x — 0, 1) can be represented as (modulo 

\x) — > {-lT\x) + \l-x) (6) 

i.e. 

Im - 1/2) ^ (-inm - 1/2) + |l/2-ni) (7) 
Following the discussion preceding equation (3), the general q-deformed state 

|0)ig|0)2g is the ground state corresponding to two non-interacting q-deformed 
oscillators In our notation a qubit state has either (a) ni = 0,n2 = 1 or 
(b) ni = l,n2 = 0. Hence 

i".-v2>.^ (S;^i:gj °). ; i-(»,-i/2)),. g];^:^-;; io), 

(8) 

So the Hadamard transformation for q-deformed state is 

Im - 1/2), (-1)"! Im - 1/2), + 1 1/2 - m), (9) 

This simphfies to: 

[F,{N,,q)a\r[F2{N,,q)al]'--^\cl>),^ 
+ [F,{N„q)a\]'-'^^[F,{N„q)alr\<P), (10) 
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where 



(11) 



ni,n2 is always or 1 so as to correspond to the qubit. It is simple to 
check that the q-number [0] is equal to the ordinary number and simlarly 
the q-number [1] equals odinary number 1. Hence the q- numbers [ni], [712] 
are always the usual numbers ni,n2- Same restrictions also apply to usual 
(i.e.undeformed) oscillators. So we restrict the hatted number operators, Ni 
and hy Ni + N2 — I where / is the identity operator. 

The Phase Shift transformation of qubit states is : |a;) — > e'^^^\x) i.e 
|0) |0) and |1) e*^|l). which in our notation is \n — |) — > e^'^^\n — |) 
where 9 is the phase shift. Now one can proceed as described in the previous 
sections. Details are in Ref.[l]. There it was shown that both the Hadamard 
and Phase Shift transformations can be realised with q-deformed qubits. 
Below it is shown that the same is possible for both the NOT gate and the 
CNOT gate. 

S.The NOT gate 

The NOT gate essentially flips a qubit, i.e. |0) |1) and |1) |0). It 
acts on a qubit as : \x) ^ \1 — x) where x — 0,1. In our notation this is 
1^1 — I) — > || — rii). For q-deformed states this means \ni — ^ || — ni)q. 
In terms of q-deformed oscillator states this becomes 

i.e. the exponents of the two creation operators are interchanged. Rewritten 
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in terms of the functions F this looks hke 

[F{N)a\nF{l - N)a\Y--\ct>), - [F{N)a\Y--[F{l - N)a\m, 

where one has used the fact that rii + n2 — 1 and followed the arguments 
after equation (136) of Ref.l, relabelled rii as n etc. Using (Ic) one gets 

[F(N)r[F(l + n- N)]'-^(a\nal)'-^\<l>), 

[F(7V)]i-"[F(2 -n- 7V)]"(al)^-"(4)"|0), (13) 

With respect to the states \(f))q, the above expression would be indistinguish- 
able from the usual " NOT" transformation if 

[F(7V)]"[F(1 +n- 7V)]^-" = [F(7V)]^-"[F(2 -n- TV)]" 

which simplifies to 

F{N) = F(l - TV) (14) 

for both n = and n = 1. 

Written in terms of its eigenvalues means 

V^i(g) _ (g"" - ng~" - ng"~^) 

This has the solution 'il)\{q) — '02(g) — '0(g) (say) for both h — Q and h — 
1. Thus the NOT gate is reahsable with deformed qubits. Moreover, the 
conditions for realisation is the same (i.e. ■0i(g) = '02(g) = '0(?)) ^ foi" the 
Hadamard and Phase Shift gates. 
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4. The CNOT gate 



The Controlled- NOT gate is a two-qubit operator where the first qubit is 
the control and the second qubit the target. The action of the CNOT gate 
is defined by the following transformations: 

|00) ^ |00) ; |01) |01) 

|10) ^ 111) ; 111) ^ |10) 

where |00) = |0)|0); |01) = |0)|1) etc. The first line of the transformation 
signifies that when the control qubit is in the "0" -state , the target qubit 
does not change after the action of the CNOT gate. The second line means 
that if the control qubit is in the "1" -state, target qubit changes value after 
the action of the CNOT gate. This may be written as (modulo constants) 
as \xy) ^ {1 — x)\xy) + x\x 1 — y) i.e. \x)\y) ^ {1 — x)\x)\y) + x\x)\l — y) 

Let the oscillators corresponding to the |a;) qubit be denoted by a, a) and 
those corresponding to the \y) qubit be b, 6^Then in terms of oscillator states 
the CNOT transformation reads: 

;|0)a- 



(N!N!)V2i^'"([fc,]![fc2]!)^/^' 



(N!N!)i/2"''"(N!N!)V2' 



n ^ ^ v"2j I i\ y^i \i\ 



where ni,n2 and ki,k2 are the eigcnvahies of the number operators corre- 
sponding to the respective oscillators with rii + n2 — l,ki + k2 — 1 and 
\(f>)a, \(f>)b denote the ground states corresponding to oscillators ai^2 and 61^2 
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respectively. Writing, 

the equation (16a) for the CNOT transformation looks like: 

\r]i)\r]2) ^ (l-ni)|77i)|772) + ni|77i)| -77) (166) 

In all subsequent discussions we shall use this form (166). However, for 
completeness, we note that the CNOT transformation (16a) can also be 
written in the alternative notation as 

11 11 11 

\ni-^a\ki--^h {'^-ni)\ni-^a\ki-- >b +ni\ni- -)a\--ki)b (166) 

i.e. 

I 2'"l 2 ~^ 2 2 ' 2 2 ~^ 2 2 

,1 I _ 1> ,h ,h . |1> ,1 ,1, ,1, 

I2 I 2 ~^ \2'"-'2'^ ' '2 2 12'^' 2 

For deformed qubits the CNOT transformation will be 

\^)q\y >q-^ (1 - >q Iv) q + " V) q (l^a) 

or in terms of deformed oscillators : 

(H!N!)V2l'^^"^([^,]![^,]!)i/2 

^ '^(N!N!)V2l'^^<^«([)ti]![A;2]!)V2l'^^^^ 
+"^(N!N!)V^'^^-(N![^i]!)V^'^^^^ ^'''^ 
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As in Ref.l, the harmonic oscillator realisations for the operators a^, and 
bq, respectively are written in terms of the two functions F and G as^^'^^: 

a\q = F{N,q)a[ ; 4, = F(l - TV, g)4 (18a) 

h{q^G{k,q)h{ ■ bl^G(l-K,q)bl (186) 

N and K are the respective mimber operators with cigcnvahics h and k and 



F{N,q) 



, G{K,q) 



Using these expressions in (176) (and relabeling rii as n and ki as k etc.) and 
suppressing the q dependence in F and G to avoid cumbersome notation one 
gets 



^ ^ ^ ^ ([n]![l-n]!)^' ^ ' ^ ^ ^ ^ ([A;]![l - A;]!)^ "^^'^ 



(l-n)F"(7V, ,)F(l-7V+n)-"-|^|^^^^ 

([nj![l — nj!)2 ([A;J![1 — /cj!) 2 

+nF"(7V, g)F(l-iV+n)^-" ("i)"^"^)'"" \<p) ^^G'^'^ {K , q)G{l-k+kf 



([n]![l-n]!)5 v v j ([1 - fc]![A;]!)^ 

(20a) 



Denoting 



(N!M!)V2'-'"« ([^illf^^lOV^'"'"^ (N!N!)V2 

F'*(7V)F(l-iV + n)^-'' = ^ 
G^{k)G{l - k + kf-^ ^ B 
G^-^{k)G{l - k + kf = B' 
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the equation (20a) becomes 

A\I3^,)B\P2,) ^ (1 - n)A\l3^,)B\l32,) + nA\l3^,)B'\ - p^g) (206) 

Multiplying both sides of (206) by {AB)~'^ gives 

\Pi,)\P2,) ^ (1 - n)\P^,)\P2,) + n\(5,,)B-^B'\ - P^,) (20c) 

Note that with respect to the states ,(20c) will be indistinguishable from the 
usual CNOT transformation (166) if B^^B' = I (the identity operator) i.e. 
B = B' or 

I (l-i^)(g-g-i) ) [ (K - I + k){q - q-^) I 
Equation (21) is an identity for both /c = and k = 1. Therefore the 
condition B — B' is always realisable in the domain of k. So the two qubit 
CNOT gate can be reahsed with q-deformed oscillators. Hence all the gates 
required for universality can also be realised with q-deformed oscillators. 
This implies that any quantum logic gate can be realised with q-deformed 
oscillators. Thus quantum computation has an alternative formalism. 

5. The possible states 

There are two possibilities as regards the arbitrary functions il'i,2, Pi,2- 

Caseil All of them are unity and hence N — N and similarly K — K. 
So (2a) just relates the opertors a, with a^, aJ.A similar argument holds 
for the operators 6, 6^ and 6g, 6j. Also from (26) we then have N = N and 
K = K. This means that at the occupation number level the deformed states 
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cannot be distinguished from the undeformed states and we are in the realm 
of standard quantum computation. Wc denote eigenvalues of the number 
operators for deformed oscillators in Case 1 hy n,k {n, k still correspond to 
undeformed oscillators); the states in Case I by D/.Then relabel rii by n etc. 
we have for Case I 

\n - 1/2) i\k - 1/2) I = \n)iai\l - n) ia^\k) ibi\l - k)ii,^ 

= |^)/ai|l - ^)/a2l^)/6i|l - ^)/62 (22a) 

where n = 0, 1; A; = 0, 1 and n — h;k — k. Note that all states have j — |.The 
possible states are: 



|00), = 


1 

' ~ 2 




1 s 

1 2) la 


1 ^ 


— 0)701 1)702 0)/6jl)/62 


(226) 


|01)7 






1 s 

1 ^/la 




0)7ai 1)702 l)/6i 0)/62 


(22c) 


|10)7 


1 

~ '2 




\\ 1 




|i)7ai|6)7a2|0)76i|i)762 


{22d) 




111)/ 


-'2 2>' 


-1^ 1 




1)701 0)7a2 1)761 0)762 


(22e) 



Case:II 



We have a general scenario if the arbitrary functions ipi{q) , (3i{q)i — 1,2 
are not all equal to unity. As these are arbitrary, let us choose ipi — ip2 — 
= (32 = f3. Then N ^ N - (1/s) In ^(g) [(26)]; and K ^ K - 
(l/s) ln(3{q). 

Hence states labelled by the occupation number are different as the eigen- 
values of the number operator of standard oscillator states and the eigenval- 
ues of the number operator of deformed oscillator states are now related by 
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n = n — (1/s) In ip{q) ; k = k — (1/s) In (3{q). This would show up in 
the Jordan-Schwinger construction. We denote eigenvalues of the number 
operators for deformed oscillators in Case II by n', k' and the states by | >//. 
So 

\n' - ll2)n\k' - 1/2),, 

= |^')//ai|l - n')na2\k')IIb^'^ - k')iib2 

= \n-{l/s)ln^l))naA'^-n+{l/s)ln%l))iia2\k-{l/s)ln(5)iii,-,\l-k+{l/s)l^^ 

(23a) 

All possible states are : 



|00),, = 1 - 


1 

2 




Iv _ 
2)//- 


1 ^ 1 




— |0)77ai 1)7/02 0)i"/6i|l)-f/62 


(236) 


|01)7, = 




1 

2 


2)//- 


1 - 2'na 




0) 77ai 1) //a2 1 ) //6i 1 0) 7/62 


(23c) 


|10)//-| 


1 
2 




Iv _ 

2)//- 


1 


Iv _ 
- 


1 1 ) //ai 0) 77a2 0) 77^^ 1) 7762 


(23ci) 


111)// 




1 

'2 


1^)11- 


l2^''"'2^ 


//6= |1) 


//ai 0) 77a2 1) 776i 0) 77^3 


(23e) 



Consistency demands the following interpretations: 

(1) For (236), ip — q^; P — q'^ i.e. the qubit state |0)77aji)77a2 corresponds 
to an oscillator occupation number n > while |0)77bjl)77b2 corresponds to 
an oscillator occupation number k > 0. 

(2) For (23c), ijj — q'^; f3 — q''~^ i.e. the qubit state 1 0) 7701 1 1)7702 corre- 
sponds to an oscillator occupation number n > while 1 1)7761 10)7762 corre- 
sponds to an oscillator occupation number k > 1. 
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(3) For {23d), tp = g""^; j3 = q'^ i.e. the qubit state |i)//ai |0)//a2 corre- 
sponds to an oscillator occupation number n > 1 while 10)7/6^1)7/62 corre- 
sponds to an oscillator occupation number A; > 0. 

(4) For (23e), il> = q^^^] j3 = q'^'^ i.e. the qubit state |i)//„i |0)//a2 cor- 
responds to an oscillator occupation number n > 1 while |l)/76i|0)/762 corre- 
sponds to an oscillator occupation number k > 1. 

Therefore we always have n > n',k > k' . ip{q),P{q) cannot be unity (i.e. 
h, k cannot be zero) because then we will have n' = h, k' = k i.e. Case I. So 
the deformed states in Case 11 can be related to harmonic oscillator states 
with occupation numbers greater than zero. 

Denote the F and G functions corresponding to the two possibilities by 

Fj, Gi and Fjj, G/7.Then 

\N(q-q-^)) \K{q-q-')) ' 

m,-r^) ) ■•''"^''''^-[ k(,-r') ) 

where we have labelled the arbitrary functions by ■1/^1,2 and /3i,2- Now the 
properties of the operators F and G have to be understood in terms of 
their eigenvalues.Then the ratio of the eigenvalues of F// and Fj is (choosing 

■01 = '02 = and /3i = = (3) 

Eigenvalue of Fn /g^>i(g) - ^'2(9) V''^ ,1. v v 
So we may write 

FH^ip'Hq)Fi (266) 
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Similarly 



Eigenvalue of Gn _ ( q^''Pi{q) - (^2{q)V''^ _ . . 
Eigenvalue of ' \ J " ^^^^ ^^^"^ 



So we may write 



Gn = (5Hq)Gi {27b) 



Thus 



' (K]![n'2]!)V2 IWa, ([fc/]![fc/]!)i/2 IWt 



Therefore 



(K]![l-n']!)V2 
([A;']![1-A;']!)V2 l'^/^^ 

Ti' 1— n' fc' 1 — fe' 1 1 

^ i/j^i/j^- P^~\n,k)i ^ i/j^p^^n^k)! (28) 



/(n, A;|n, A;)j 

= ^5(g)/5^(g) (29) 

So the right hand side of (29) is a function of q only. For ip{q) = 1 and 
P{q) = 1, one cannot distinguish between the two cases at the level of the 
scalar products between states. However, if the arbitrary functions are not 
unity then these scalar products are distinct from each other and this might 
be useful at the level of experimental realisations or consequences. Moreover, 
note that one can choose the arbitrary function ip{q) for the control qubit 
to be the same as that in the Hadamard, Phase shift and NOT gates. As 
the CNOT gate is a two qubit gate, a different function I3{q) is taken for the 
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target qubit. So with two arbitrary functions all the three gates required for 
universality can be constructed with q-deformed qubits 

6. Conclusion 

Therefore quantum computation admits of an alternative formalism where 
q-deformed oscillator states can be used to construct qubits. This has been 
estabhshed here with the realisation of the CNOT quantum logic gate with 
q-deformed oscillators. Thus this realisation is possible for all the quantum 

logic gates required for universality. Hence all quantum logic gates can be 
realised with q-deformed qubits. The existence of additional parameters will 
enable comparison between different experimental scenarios using the usual 
scheme and the alternate scheme. This requires further study. 

Another aspect where the new formalism might prove useful is the realm 
of quantum error correction. The conventional way to correct computer error 
is to use redundancy. More than one clement is used to denote the same 
bit. For example, consider two atoms A and B and use the doublet AB to 
store the same bit (of information). Thus one has the state |00) = |0)|0) (or 
the state |11) = |l)|l)).An error changes the state of only one atom, i.e. one 
now has any one of the states |01), |10) instead of |00) (or the states |10), 
|01) instead of |11)). For type II states, |)//, we have n' — n — Hnip. It is 
simple to check that ip = q"' for n' = and = q'^~^ for n' = 1. Therefore 
occurrence of an error is reflected in the change of the arbitrary function and 
consequently in the matrix elements n{n'\n')ii. So the arbitrary functions 
provide an additional leverage and can possibly be used in detecting and 
regulating errors. This also requires further investigation. 
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